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ABSTRACT 


This paper proposes the adaptive particle swarm optimization (APSO) 
technique to control the active and reactive power produced by a variable 
wind energy conversion system and the exchanged power between the electric 
grid and the system during a voltage dip (VD). Besides, to get the variable 
speed wind energy maximum power, a maximum power point (MPP) 
methodology is utilized. The system under study is a 5 MW wind turbine 
connected via a gearbox to a doubly-fed induction generator (DFIG). The 
DFIG stator is branched directly to the electrical network, while the Back-to- 
Back converters couple the rotor to the grid. The decoupled vector control of 
the rotor side converter and the grid side converter is established primarily by 
a conventional proportional-integral (PI) and a second level by an intelligent 
PI whose gains are tuned using the proposed control. The performances and 
results obtained by APSO tuned PI controllers are analyzed and compared 
with those attained by classical PI controllers through the 
MATLAB/Simulink. The superiority of the advised technique is examined 


during a two-phase short-circuit fault condition and confirmed by the reduced 
oscillations. 
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1. INTRODUCTION 

To have environment-friendly energy, the industry is perpetually looking for new techniques and methods 
whose purpose is to produce electricity from natural resources (wind, solar rays). Wind energy is generated by 
converting the kinetic wind energy. Recently, this energy has known great importance in the countries which invest in 
the sector, thanks to the immense progress made in scientific research. The most common configuration of a variable 
speed wind power plant (WPP) is one with a doubly-fed induction generator (DFIG), whose stator is linked directly to 
the grid, as indicated in [1]. However, the back-to-back converters connect the rotor to the grid, allowing for 
bidirectional power transfer, as seen in Figure 1. The fundamental benefit of this architecture is that the converters 
allow for the transmission of a percentage of the total system power, which reduces losses in power electronic 
components [1]. The generator control is based on the stator or rotor flux oriented (SFO or RFO) vector control 
approach with a conventional PI controller. Laafou, et al. [2], the rotor side converter (RSC) regulates the stator active 
and reactive power generated by the DFIG by controlling the rotational speed and the rotor currents using 
conventional PI controllers. However, the grid side converter (GSC) control maintains the DC link voltage constant. 
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Figure 1. Synoptic scheme of the wind energy conversion system (WECS) connected to the grid 


In Ahyaten and Bahaoui [3] have developed and simulated the wind power plant using a wound 
rotor induction generator DFIG connected via back-to-back converters GSC and RSC to the non-ideal 
electrical grid. Moreover, the authors carry out the control of the pitch angle. Likewise, [4] has studied 
different control schemes for the RSC which is utilized in the wound rotor induction generators. Examples of 
such control schemes are flux oriented (FOC), voltage oriented (VOC), direct torque (DTC), and direct 
power control (DPC). Additionally, authors in [5] have introduced the modeling of the components of DFIG 
based on a wind energy conversion system. The back-to-back converters and a DC link capacitor connect the 
rotor to the electrical network. The suggested technique is based on the feedback controller to pilot the 
dynamic performance of the DFIG and control the maximum produced power. Elazzaoui [6] has modeled 
and simulated the wind energy conversion system using a doubly-fed induction generator of 3 MW. The 
control scheme is based on conventional PI controllers where the gains are set by the pole compensation 
method. The previously mentioned articles show that the researchers did a great work modeling the wind 
energy conversion system and established the different control schemes, principally based on the classical PI 
controller. As investigated in [7], this type of control has limitations as its performances depend on 
calculating the conventional PI controller where the gains are constant. As well, it is still suffering from 
dependence on the parameters of the DFIG machine. These gains can probably change following an increase 
in temperature during the generator operation or when a defect occurs in the generator. The voltage dip (also 
denoted voltage sag) is a grid perturbation and defined as a short-duration fall in the RMS value of the 
voltage induced by a short circuit fault, which influences the PI controller performance [8]. Rafiee, Rafiee, 
and Aghamohammadi [9] has studied the effects of the voltage dip that occurred in the DFIG stator and rotor. 
Besides, the DFIG is extremely susceptible to grid disturbances such as voltage sag owing to short-circuit, as 
indicated in [10], [11]. This fault can damage the converters and provoke the wind system disconnection 
from the grid [12]. The authors have studied the low-voltage ride-through (LVRT) methods, which allow 
system protection. However, in [13], authors have proposed a comparison of the PI controller and the sliding 
mode technique for the grid side converter under voltage dip. To deal with the drawbacks of the traditional PI 
regulator, the researchers have proposed various intelligent controllers. One of the most extensively used 
artificial intelligence methods is particle swarm optimization (PSO). PSO is applied to get optimal gains of 
controllers. Taleb, Cherkaoui, and Mahfoudh [14] have applied particle swarm optimization (APSO) to 
enhance the performances of the PI controllers. The latter is utilized to control the rotor side converter. A 
comparison of the suggested controller to the conventional PI controller is established and offered. Laina et 
al. [15] have proposed a comparison analysis of PSO with sliding mode control (SMC). The authors compare 
PSO-PISMC, PSO-ISMC, and the traditional PSOSMC to show that the PSO-PISMC presents a faster 
transient response and a lower tracking error for a wind turbine system utilizing the DFIG. Likewise, the 
authors [16] have examined the overshoot anomaly of the stator and rotor currents when there is a 
perturbation in the grid and have proposed the PSO algorithm to tune the parameters of PI controllers for 
reducing the currents peak value to protect the back to back converters. In this research, we propose the 
APSO for tuning the gains of the PI controllers, improving the response time of the system, reducing the 
overshoot, and ensuring better tracking of the reference under a varying wind velocity and the grid voltage 
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dip produced by a two-phase-to-ground fault. Thus, the major contributions of this article are summarized is 

being as: 

— Modeling and simulating the proposed optimization technique for the WECS, which debiting in the grid 
under-voltage dip using MATLAB/Simulink. 

— Results are compared to those found in the WECS using conventional controllers PI. 

— The suggested method performance is analyzed in terms of the convergence of powers and the DC bus 
voltage. 

The following is how the rest of the article is structured: The study of mathematical models of the 
wind power plant is covered in section 2. The study of vector control using a traditional PI for the wind chain 
is the focus of section 3. Section 4 highlights the APSO. Section 5 tackles the results and comparison 
between the two techniques. 


2. MODELING OF THE WIND ENERGY SYSTEM 
2.1. Wind turbine modeling 
The aerodynamic power produced by the turbine is expressed as shown in (1), [17]: 
Pry = — Cp A R2. Vw? l 
Tu = 5: Cp (A, B).p. T.R. Vw (1) 
Where Cp is the turbine efficiency performance coefficient of power. This coefficient can be estimated using 


(2). It depends on the ratio between the linear speed of the blades and the wind speed (A), which is expressed 
as shown in (3), and the pitch angle 6 [18], [19]: 


Cp(B, A) = [0.5 — 0.0167. (B — 2)]. sin ee) — 0.00184. (À — 3). (B — 2) (2) 
R. Ory 
j= Row (3) 


Supposing that the overall mechanical dynamics of the system are brought back to the turbine shaft, the 
following equations describe the model [6]: 








Jtot- dt E Qmec = Tg- Tem (4) 
= TTu n Omec 
T= F and Gp= o (5) 


Where Jtot is the overall inertia of WECS, Tru is the turbine torque, T, is the gearbox torque, Tem is the 
electromagnetic torque of the DFIG, and fis the overall viscosity coefficient of friction. 

The MPPT algorithm is utilized to obtain the maximum power produced by the turbine. Its principle 
consists of controlling the electromechanical torque of the DFIG to adjust the rotation speed by maximizing 
the power. In this study, the MPPT with the speed regulation is selected based on a typical PI controller. The 
pole compensation approach is employed to determine the gains of the controller. The following are the 
expressions for these parameters [20]: 


1 . l . 

Kimppt = m With T is the system time constant (6) 
—Kimppt Jtot 

Kpmppt = E (7) 


2.1. Modeling of the doubly-fed induction generator 
The following electrical equations can be used to define the dynamical model of the examined DFIG 
using the Park transformation [3], [21]: 
= R.,..i WWstad _ 
i ENS dt = vee Ystad = Leta. Ista_d + Ly. lrot_d 
Wsta_q = Leta- Ista_q + Lin. lrot_q 
Prota = Lrot-lrot_a + Lm: Ista_a 


. dWsta_q 
Vsq = Rsta Ista.g + di + Ws. Wstad 
oy . dWrotd 
Vrot_d z Rrot: lrot_d + dt — Wrot: Wrot_q 


; dŲrotq 
Vrot_q = Rrot- lrot_q + + Wrot: Prota 
dt 


(8) 


rot_q — Lrot: lrot_q + Lm: Ista_q 
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The stator and rotor voltages are V, and Vio , the stator and rotor currents are ista and irot, Ysta and 
Wrot are the stator and rotor flux linkages, Rsta and Riot are the stator and rotor resistances. Ws and Wrot are the 
stator and rotor angular frequencies. The stator, rotor and magnetizing inductances are Leta, Lrot and Lm, 
correspondingly. The electromagnetic torque can be given as [21]: 


t= "p. - Crot q: Ystad = teed: Wsta.q) ; With p is the number of the DFIG pole pairs (9) 


3. VECTOR CONTROL APPLIED TO DFIG-BASED WIND ENERGY 
3.1. Application of the indirect vector to control the rotor converter side 

The orientation of the stator flux vector, which is oriented along the d-axis, is used to control the 
RSC. The phase-locked loop (PLL) strategy is employed to calculate the electrical network voltage (Vs), and 
the resistance of the stator windings is ignored. 


Wsta_d = Wsta ’ Wstaq = 0 (10) 
Vsa = 9, Vsq = Vs = Ws Wy, (11) 


The rotor currents can be expressed as shown in [6], [21]: 





lrot_d = [Viota + 8: Ws (Tiss a m”) lrot_q | /[Rro T (Liot — : s| (12) 
lrot_q = Veet — 8 Ws (L E m") lrot.d — ae 1/[Rrot + (Liot -_ m”) . s| (13) 


The electromagnetic torque can be rewritten as: 
_ Liig 
Tem = —Pp. L -lrot_q: Ysta (14) 
sta 


The goal is to control the stator power independently pumped into the grid. The following equations 
can be used to express the stator active and reactive power [6], [21]: 


Lm Lm 








P; = —Vsq: D, irot_q = -Vs 7, Irot_q (15) 
Véq Lin . VE Lm . 
— = — ] = — V. — i (16) 
Qs Ws-Lsta sq Lsta rot_d Ws-Lsta 3 Lsta rot_d 


An open loop regulates the powers. The rotor currents (irot g> lrot_a) are, nevertheless, controlled by 
a closed loop. Unlike, the current i; q is obtained from input power reference to control active power. The 
reactive power is controlled by regulating it, g, which is computed from the Q,’. These current references 
can be written as in [15]: i 








. = Lsta 

irot_q = ts lin We . Pe (17) 
“* feat a Lsta xo V: 18 
lrot_d ~~ Fa Ve . ( Qs Ws-Lsta ( ) 


The voltage references are expressed as shown in [22]: 
* .* . 1 1 
Vrotq = Trot ee [Kp-rsc1 + Ki-rsc1: 3] + era t Vs (19) 


* -x . 1 
Vrot.d = liia E lavd: [Kp-rsc2 a5 Ki-rsc2:- 7] F erq (20) 


Where 
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_ ae = ha N Bethy => <1 Vela 
Crd — 8. Ws. Ca — =.) ‘lrot_.d) €rg = 8. Ws. (Dei — ae -lrot_q , Vs = 8- (ees (21) 


Considering T;sc the system time constant for the RSC controller, the PI gains can be determined 
applying the pole compensation technique as shown in: 


—_1 Da _ Kprsc-Rrot 

Kp-rsc1,2 T an (Lrot = L ) And Ki-rsc1,2 = Len2 (22) 
rsc sta (L _2m ) 
rot Lsta 


3.2. Control of the grid side converter 

The GSC is responsible for regulating DC-link voltage as well as reactive power exchanged with the 
grid. The grid voltage is oriented belong the q-axis and is determined by the PLL technique to elaborate the 
GSC control strategy. The grid voltages can be expressed as: 


Vea = 0 And Vyq = V; (23) 


The voltages and powers can be simplified as shown in [6], [23]: 


Vasca = —[Re + Le s]. igsca + We: Lt- igseg (24) 
Veseq = — [Ret bes Fiza = Orbi lasca F Vg (25) 
Ps. = Ves lesoq (26) 
Qg = Vg: Igsca (27) 


The relation between the powers of converters can be expressed as shown in [6]: 

Vac: lc=Pesc — Prse (28) 
The following is how the GSC power can be deducted: 

Pe = Pesce = Vac-lc + Prse (29) 
Where Prsc is the rotor side converter power, which can be defined as: 

Pesc = Voc: Irsc (30) 
So, the DC-link power (Pje) can be expressed as: 

Pi. =Vpc. it eo) 
From (26)-(29) the references of the grid currents can be derived as shown in [6]: 


= ok 1 * ax shiek = Qs 
Isscq = Ve (Voc-ic + Prsc) ; loscd = v (32) 


It is assumed that the reference of reactive power is null to have a unit power factor. So, we impose a null 


* 


direct current grid reference (igs¢q = 0). The voltage references are expressed as shown in [23]: 


Via = [ižsca — igscal-[Kp-esce + Ki-gsc2- <] + Egscq (33) 

Via = [iġsca — igscąl: [Kp-gsc1 + Ki-gsc1-5] — egsca + Vg (34) 
Where, 

egscqg = Og- Lf- igscq) Egscd = Wg: Le-Igsca (35) 
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The expressions of the grid currents can be deduced as shown in [24]: 


. 1 x n 
. S 1 x 
loscd = ee Les] ; (Veg T We. Le. loscdq m Voc: Sq ) (37) 


Where Sa and Są are the switching states computed by Park transformation. To keep the DC-link voltage at 
its reference, the PIpc controller is utilized. As a result, the PI controller parameters are as shown in: 


Kp-pc = 2.§ w.c and Ki-pc = w*.c ; Where & is the damping coefficient (38) 


The currents igscq and igsca, flowing through the RL filter, are regulated by a PIcsc controller used in the inner 
loop. Considering the time constant Tgsc of the controlled system, the PI controller gains are: 


Lf 


Rf 
K,_ =— ; K_ = 
p-gsc1,2 Tec » **1-—gsc1,2 


ae 39 
Teso (39) 


4. ADAPTIVE PSO BASED PARAMETERS SELF-TUNING 
4.1. Overview of PSO 

PSO is an evolutionary meta-heuristic computation methodology inspired by examining the 
behavior of bird flocking or fish schooling when they are hunting for food or keeping away from opponents 
[15], [16]. Kennedy and Eberhart devised it for the first time in 1995. The PSO is an artificial intelligence 
methodology based on decentralized and self-organized collective behavior. The population candidate 
solutions are used in the PSO methodology to find an optimal solution for a problem [18]. A fitness function 
determines the optimality level. The PSO method works on the notion that each particle is looking for the 
objective in D-dimensional space. The other particles will immediately adopt the best and ideal position that 
has been discovered [25]. 


4.2. PSO mathematical model 
Each particle (xij), randomly defined by initial velocity and position, tries to explore the research 
space through iterative test positions to achieve their objective following the expressions given by [26], [27]: 


Vij (t =F 1) = W. Vij (t) F C4. rı. (Be = Xij ©) F Cj: r2. (Goest = Xij ©) (40) 


Where Vjj(t) represents the particle im velocity with a j dimension at iteration t, x;j(t) represents the particle im 
position with a j dimension at iteration t, Pres is the best prior position of im particle, Gvest 1s the best particle 
among all the population, W is the inertia weight factor, (Cı and C2) are the acceleration constants, (rı and rz) 
are the random integers between [0-1], and n is the swarm size. 


4.3. Tuning PI parameters using APSO 

The PSO algorithm calculates the optimum PI parameters [Kp, Ki]. Each particle swarm targets a 
minimum of the fitness function to guarantee the system stability in a closed-loop. The integral absolute error 
criterion (IAE) is employed as a performance index to minimize the fitness function in this research. The 
following expression defines the IAE [28], [29]: 


IAE = f, |e(t)|. dt (42) 


The signal error e (t) is computed by a difference between the reference input and the actual output 
value [27]. The weighted sum of errors, as given in (43), is determines the objective function. The target 
function in (44) is associated with the IAE performance criterion. By minimizing the fitness function 
employing the adaptive weighted PSO, the optimal parameters are obtained with fast convergence of the 
algorithm [28]: 
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FD=DE, w len! (43) 
Obj(t) = A F(t) dt (44) 


Where d is the dimension of the problem. The weighted factors (wi) are positive and related to the fitness 
function, which their sum is [29], [30]: 


Ziwi = 1 Withi= 1...5 (45) 


4.4. Design of the algorithm 
The steps of the searching procedure for the PSO technique are listed as shown in [25], [27]: 

— Step 1: Randomly generate an initial population. 

— Step 2: Once a certain number of iterations have been completed, the algorithm is terminated. 

— Step 3: Evaluate each particle's objective function and record each particle's best prior position (Pi) as 

well as the global best position (Pg). 

— Step 4: For each particle, update the improved velocity of formula (40) and the position of formula (41). 

Verify the velocity constraint conditions by using the following expressions: 


Vnaz if Vi (t + 1) > Vraz 
Vmin if Vi (t + 1) < Vmin 


Verify the position constraint conditions of particles as shown in: 


Xmax if x(t T 1) > Xmax 
x(t + 1) = x(t + 1) if Xmin < V(t + 1) < Xmax (47) 
Xmin if Xi (t + 1) < Xmin 


— Step 5: Go back to step 2. 


5. DISCUSSION OF SIMULATION RESULTS 

Matlab/Simulink software is employed to model, simulate, and perform the wind energy conversion 
system. Figure 2 presents the wind profile. The DC bus voltage reference is constant and is equal to 1200V 
(Vpc*=1200 V). Figure 3 shows the rated mechanical power reference, which is obtained by using the MPPT 
algorithm. It can be seen that when the wind profile equal to 12.5 m st, the rated produced power is 5 MW. 
The doubly-fed inductor generator is driven by the speed rotation of the multiplier output. The mechanical 
speed, shown in Figure 4, is obtained by the MPPT algorithm. The reactive power references Q,* and Q,* are 
set to O (VAR) to ensure a unity power factor. The pole compensation method is applied to calculate the 
parameters of PI controllers. These parameters, given in Table 1, are tuned by using the APSO. The 
configurations illustrated in Table 2 are chosen to set up the PSO program. The obtained results are compared 
between the conventional PI and the tuning APSO controller PI. In this paper, we analyze the active and 
reactive power, which are generated by the (DFIG) during the voltage dip. Figure 5 shows the grid voltage 
signals during a voltage dip, which occurs at 0.3 (s) and holds out 200 (ms). However, Figure 6 shows the 
quadrature grid current computed by the PLL. As you can see, the voltage drops by 80 % due to a short 
circuit between two phases. Figure 7 illustrates the pattern of the reactive power during the voltage dip. The 
response time is ameliorated, and the ripples are significantly reduced by utilizing APSO. It can be observed 
by comparing the results shown in Figure 8 and Figure 9, that the active power produced by the DFIG (Ps- 
PSO) perfectly follows the reference (Ps) than (Ps-PI). The response time is considerably improved to 4.3 
(ms) from 10 (ms), and the static error is null in steady-state. When a grid fault happens at 0.3 (s), the PSO 
reduces the ripples with an ameliorated time response. The DC-link voltage computed by the traditional PI 
and the smart PI are shown in Figure 10. The APSO approach reduces DC link voltage overshoot to 1399 V 
from 3452 V and increases response time, as seen in Figure 11. The performance comparison is summarized 
and presented in Table 3. The comparative findings of the two proposed controls are shown in this section. 
Table 4 shows the characteristics of the wind power plant components used in the simulation. 
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Figure 5. Grid voltage Vg (a, b, c) during the VD 








Figure 6. Quadrature grid voltage during the VD Figure 7. Reactive power during the voltage dip 
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Figure 8. Active power during the VD Figure 9. Active power during the VD-zoom 





Figure 10. DC link voltage during the voltage dip Figure 11. DC link voltage during the VD-zoom 
Table 1. Gains values of PI controllers Table 2. Parameters of PSO code 
Paramaeters PI without APSO PO with APSO Parameters Value 
Kp dc 1,848 1,954 Population size 20 
Ki dc 396 322,226 Number of parameters 10 
Kp_rsc1/2 1,444 0,672/0,1067 N° of iterations 20 
Ki-rsc1/2 0,2369 0,0657/0,85 W 0,9 
Kp_gsc1/2 200 176,288/123,929 C=C, 2 
Kigsc1/2 5e+4 4,71e4°7/ 4,1e+” 
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Table 3. Performance values 
Performance Without APSO With APSO 


Active Power Response time 10 ms 4,3 ms 
Static error 0,242.10° 0 
Peak overshoot 10.10° 7,79.10° 

DC link Voltage Peak overshoot 3452 V 1399 V 


Table 4. Set of parameters used in the simulation 


Components of WECS Parameters Symbol Value 
Turbine Radius of blade R 51,583 m 
Coefficient of multiplier Gp 47,23 
Total moment of inertia Stot 1000 kg.m? 
DFIG DFIG rated power Ps 5 MW 
Stator leakage inductance Lista 1,2721 mH 
Rotor resistance Riot 1,446 mQ 
Rotor leakage inductance lot 1,1194 mH 
Mutual inductance ee 0,55187 mH 
Stator line to line voltage Vs 950 V 
Capacity DC-link capacitance C 4400 uF 
Filter RL Resistor of the filter Re 20 Q 
Inductance of the filter Le 0,08 H 


6. CONCLUSION 

The WECS is modeled and simulated in this research using changing wind speed and grid voltage 
drop situations to verify if the recommended method is superior. The DFIG stator is directly coupled to the 
power system, and the rotor is connected to the grid via back-to-back converters. To modify the gains of PI 
controllers, intelligent PI controllers based on the APSO algorithm are created. The results of the two 
techniques are compared, and a 5 MW doubly-fed induction generator is utilized to test wind chain energy. It 
is observed from the simulation results that the APSO tuning technique not only provides satisfactory and 
enthralling results but adjusts the PI controller gains to the studied system under severe conditions, notably 
when the wind speed changes suddenly or when the system settings change. When compared to a traditional 
PI controller, the APSO approach provides superior performance in terms of improved response time and 
lower static error. Besides, it is a good technique for tracking the wind profile variation, which is considered 
stochastic and variable energy. Furthermore, the APSO reduces oscillations induced by the voltage dip in 
RMS grid voltage substantially better than indirect vector control. 
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